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1 Introduction 



The purpose of this study is to examine the functional dependence of sea 
ice brightness temperature on ice thickness at frequencies measured by the 
Advanced Microwave Scanning Radiometer on EOS (AMSR-E). That is, to 
generate a set of brightness temperature-thickness (T&-/t) curves. Several 
studies have demonstrated a relation ship between ice thickness and SSM/I 
or AMSR-E brightness temperatures ({Martin et all 120041 ; iNaoki et all . 120081 ; 



Kwok et all 120071 : iHwang et all 120071 ) with at least o ne attempt to use th is 



relationship for the purpose of ice thickness retrieval (jMartin et all . [2004) 



This relationship, however, is not a direct one in the sense that, all 
other things begin equal, changes in ice thickness will not necessarily affect 
microwave emissions at AMSR-E frequencies. Even at 6.8 GHz, the pen- 
etration depth is far smaller than all but the thinnest of ice sheets. The 
relationship is caused by the fact that thinner, newer ice, has different phys- 
ical properties than those of older, thicker ice. In particular it tends to be 
more saline, especially in the top layers. 

Three processes en sure that newer, thinner ice i s, on average, more saline 
than older, thicker ice dTucker et all Il992l: lEickenl . Il992l: IWeeks and Acklevl . 
1985 : Vancoppenolle et al. . 20071 : Ehn et al. . 2007 ). especially in the top lay- 



ers. As new ice is formed, most of the salt gets expelled, except for a small 
amount that is included as pockets of highly saline brine. The faster the 
ice grows, the more brine is included. Since thin ice conducts heat more 
quickly, it will grow faster than thicker ice on average. As a growing ice 
sheet cools, some of the water in the brine will freeze. Since ice is less dense 
than brine, increasing pressure in the brine pockets will cause some of the 
brine to be expelled from both the top and bottom of the ice sheet, produc- 
ing the characteristic 'C'-shaped profile of first-year ice. Finally, as the ice 
ages, thawing will produce channels through which much of the brine can 
drain. 

The causal relationship between salinity and brightness temperature is 
similarly circuitous, although it is a more definite one. Higher salinity pro- 
duces a higher effective permittivity whose main effect, for optically thick 
ice, is to increase the polarization difference or second Stokes parameter. It 
is not just the salinity that affects the microwave signal. Changes in tem- 
perature will have a similar effect as well as a similar cause. That is, warmer 
ice will tend to have a higher effective permittivity and thinner ice will tend 
to be warmer, again because of heat conduction. 

We use parameterised empirical relationships between ice thickness and 
salinity. Model profiles are fed to ice emissivity models to determine ideal- 
ized relationships between brightness temperature and ice thickness. 
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Figure 1: Data flow diagram describing the model set up, namely how the 
other parameters are derived from ice thickness. Rectangles are processes, 
ellipses are data. The salinity-thickness relation and the salinity model are 
described in Section T2.2I The thermodynamic model is described in Section 



2 Model set up 

We use a layered, plane-parallel, radiative transfer (RT) model called Mi 



crowave Emissivity Model for Layered Snowpack (MEMLS) (jWiesmann and Maetzlerl . 



1999al lbh. The emissivity calculation can be divided into three main modules 



First, effective relative permittivities must be calculated from the physical 
properties of each layer-chiefly, temperature and salinity. Second, scattering 
coefficients must also be calculated from values of correlation length. Finally, 
these permittivities and scattering coefficients are simultaneously integrated 
through the depth of the ice sheet. Before we can feed these quantities into 
the model, however, we must take account of the numerous and complex 
interrelationships between the different inputs, particularly since we are in- 
terested only in the functional dependence of brightness temperature on a 
single bulk property, i.e. ice thickness. For example, because of changes 
in relative brine volume at different temperatures and salinities, correlation 
lengths of the brine inclusions will be affected by these other two properties. 

The emissivity model is described in the data flow diagram shown in 
Figure [5] including data inputs and parameters. The basic modules have a 
dashed outline. As can be seen, it is a quite a complex model comprising 
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Final result 



Figure 2: Detailed data flow diagram describing the complete ice emissivity 
model. Rectangles are processes, ellipses are data. The "mixture model" is 
describe d in Section 12.11 Ot h er com pon ents are described elsewh e re, see for 



instance 



Mills and Hevgster ( 201ll ) and Wiesmann and Maetzler ( 1999a! ) 



The boxes labelled, respectively, "SnelPs law", "Fresnel equations", "Inte 
grate layers" and the box with the equation for the tra nsmission coemcent 
corres pond to equations (3), (1-2), (5-6) and (7- 

(201l|). The equations for fresh- water ice permittivity are contained in 



Huffordl (1199111 ■ For other eq uations in the permittivity module, a good 
reference is llJlabv et al.1 (jl98fil l. 
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many separate yet interacting components. The model set up is shown in 
Figure [TJ 

Both ice thickness and the salinity profile (function of salinity with 
depth) will be affected by past weather conditions. In addition, the cur- 
rent temperature profile will be determined by both the ice thickness and 
prevailing weather conditions. For this reason, we have designed a simple 
thermodynamic model which can be used both to determine the temperature 
profile, and in a crude fashion, model ice growth. This model is described 
below, along with the other components. 

The calculation of effective permittivities is described b elow. For scatter- 

ing coe fficients, we use the empirical model within MEMLS (IWiesmann and Maetzlerl . 
1999b). We relate the correlation length to brine volume as follows: 



l = k. 



(1) 



where Iq is a constant and Vb is relative brine volume. For the case of ran- 
domly positioned, oriented cylinders, the correlation length in the transverse 
direction would be related to brine volume to the first power. Fo r the case 
of ra ndomly positioned spheres, it would be to the third power. (jMaetzlerl . 
19971 ) We take the second power, thus the brine inclusions are somewhere 
between a cylinder and a sphere. 

The radiative transfer model will not be described in this report since 
it and similar models are covere d in o ther publications — see fo r instance, 
Wiesmann and Maetzler ( 1999a! ) and Mills and Hevgster ( 201ll ). We will 
note, however, that there is some descrepancy in the two models contained in 
the aforementioned publications as to when and how (by simply eliminating 
it or taking the absolute value) to drop the imaginary component of the 
reflection coefficients. In MEMLS, for instance, only the real part of the 
complex permittivi t y is u sed in the calculation of reflection coefficients. In 
Mills and Hevgsterl (|201ll ). by contrast, the imaginary part is carried right to 
the end of the calculation and real reflection coefficients produced by taking 
the absolute value. Using one method or the other can mean a difference of 
as much as 10 K. We have modified the MEMLS code so that calculations 



of refl ection coefficients match the model described in iMills and Hevgster 
(mil). When scattering is neglected in MEMLS, the two models agree to 
within half a Kelvin. 

The viewing angle is set at 55 degrees, or the same as AMSR-E. 



2.1 Model for effective permittivities 

The calculation of effective relative permittivity is one of the most important 
steps in the model. It is also the one of the least understood and most 
fraught with uncertainty. To calculate the effective permittivit y, we use a 
semi-e mpirical model combining the empirical mixture mo d el of IVant et al. 



( 19781 ) and the theoretical one from Sihvola and au Kong ( 1986 ) based on 
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a. Real permittivity at 0.1 GHz b. Imaginary permittivity at 0.1 GHz 




c. Real permittivity at 4. GHz d. Imaginary permittivity at 4. GHz 

6.0 f 1 J 0.20 [' ' ' / 




3.0 3.5 4.0 4.5 5.0 5.5 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

e' (Vant et al. 1978) e" (Vant et al. 1978) 



Figure 3: Comparison of the mixture model of Sihvola and Kong based on 
de-polarizability and the empirical mixture model of Vant using random 
values of temperature in the range of 251 to 271 K and salinity in the range 
of to 15 psu. 



Table 1: Correlation between Vant and Sihvola mixture models of effective 
permittivity for both real and imaginary components. 



V 


r(e') 


r(e") 


0.1 


0.994 


0.929 


0.2 


0.994 


0.928 


0.4 


0.995 


0.942 


0.8 


0.995 


0.974 


1.0 


0.996 


0.981 


2.0 


0.995 


0.984 


4.0 


0.995 


0.983 
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Figure 4: Correlations with wavelength of regression coefficients between 
empirical and theoretical mixture models. 
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Table 2: Fitted coefficients for semi-empirical mixture model. 



Coefficient 


value 


Cl 


0.154 


C2 


1.005 


C3 


-0.464 


C 4 


-0.036 


C5 


11.011 


C6 


1.409 



the low- frequency limit. Vant et al. give an empirical mixture model based 
on fits of measured real and imaginary permittivity to brine volume: 

e* = aV b + b (2) 

where e* is the effective permittivity (real or imaginary), Vf, is the rel- 
ative brine volume and a and b are constants. Different constants have 
been acquired for frequencies between 0.1 and 4 GHz. To extend the fre- 
quency range of the Vant models, we combine it with the following theo- 
retical, mixture model for or iented needles based on the low-frequency limit 
( Sihvola and au Kong . 19881 ): 



k V6(e2-ei)ei/[€i + P(e 2 -ei)] 
e *// " 61 + i _ PVb{e2 _ ei ) [ £l + P(e2 _ ei)] w 

where e\ is the permittivity of the host material (ice), e 2 is the permittivity 
of the inclusion material (brine), and P = 0.5 is the de-polarisation factor. 
It was observed that the results from ^ and ([3D are closely correlated as 
shown in Figure El however they differ by both a constant factor (in the 
real part) and a coefficient (in both the real and imaginary parts). These 
parameters were found to vary closely with the wavelength, as shown in 
Figures [H 

Thus, the combined model is given as follows: 

tsemi = (ci/V + c 2 )e' eff + c 3 /u + c 4 + (c 5 /v + c 6 )e" ff i (4) 

where e e ff = e' e jj + ^ s ^ e e ff ec ti ve permittivity from Equation ([3]), v 
is frequency and c\ through cq are fitted constants. Note that there is no 
constant term for the imaginary part, forcing the value to zero for zero brine 
volumes (pure ice is almost a perfect dielectric). We compare the derived 
semi-emprical model with the Vant model for a range of frequencies between 
0.1 and 4. GHz, a range of salinities between and 15 psu and a range 
temperatures between 251 and 271 K in Figure [5j The fitted coefficients, c\ 
through C6 are given in Table [2j 
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Figure 6: Different salinity-thickness relationships compared. Shading are 
residuals for the Cox and Weeks models. 



Finally, a note on brine volume. In Mills and Hevgster ( 201ll ). brine 
volume was given as: 

Vb = 4~ (5) 

where Sb is the brine salinity which, because of freezing point depression, can 
be calculated from temperature only. This is a rather crude approximation. 
A more accurate formula is as follows: 

V b = -= S » - _ (6) 

bbPb - zpb + bpi 

where pi and pb are the densities of ice and brine, respectively. This formula 
is found by writing the ratio of brine volume to ice volume in terms of density 
and salinity and solving. It directly connects the calculation of brine salinity 
(an accurate function of temperature only) with relative brine volume. 

2.2 Salinity model 

Figure [6] shows plots of different salinity-thickness (S-h) cur yes which have 



been d erived empirically and from a simple model (see below) . I Cox and Weeks 



provide the following piece- wise, linear fit of bulk salinity to ice thick 
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5 10 15 20 

salinity 



Figure 7: Model salinity profiles for different bulk salinities in 0.5 psu incre- 
ments, starting at 0.5psu and ending at lOpsu. Ice depth is normalized. 
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ness (solid curve): 



S n . 



14.24 - 19.39/i; h < 0.4m 
7.88 - 1.59/i: /i > 0.4m 



(7) 



where S ave is bulk salinity and h is ice thickness in metres. This is the only 
S-h curve that will be used in the study. Rath e r tha n trying the emissivity 
models with different curves, Cox and Weeks! (1974) provide residuals for 
the curves. We will use these residuals, along with those supplied with the 
Vant effective permittivity models to compute an error bound for the final 
results (see below). 

The other five curves are as follows. Naoki et al. (2008) generated an 
exponential fit (short-dashed curve) of surface salinity (top 10 cm) based on 
measurements conducted in the Sea of Okhotsk: 



S = exp(m/i + k) 



(8) 



where S is th e salin ity, h is ice thickness and m and k are constants. 
Hevgster et al. (120091 ) ge nerated an ex ponential fit of ice core samples taken 
from the Weddell Sea (jEickenl . [l992j). Coefficients for equation (|HJ) are 
m = —3.55 and k = 3.06358 for the Naoki model and m = —0.275 and 
k = 8.04 for iHevgster et al. (2009). Notice that the Naoki curve roughly 
matches the left section of the Cox and Weeks curve, while the Mills and 
Heygster curve roughly matches the right section. 

The dotted curves show results from the thermodynamic ice growth 
model (described in the next section) for two weather scenarios: fall freeze- 
up and a winter cold snap. These provide a reasonable range for thin ice 
(below 30 cm), however they over-estimate salinity for thick ice. 

Two trials will be performed: one in which the ice is represented as a 
single layer, and one with a parameterized salinity profile. For the param- 
eterized salinity profile, we use the 'S'-shaped model from lEickenl (j 19921 ). 
but with the caveat that at value s for bulk salinity belo w 2psu, it reverts 
to a flat profile in accordance with iGranskog et al. I (|2006h . This is shown in 
Figure [7) 



2.3 Thermodynamic model 

As mentioned in the introduction, it is not just ice salinity that lowers 
the brightness temperature of thin ice, it is also the temperature: thin ice 
conducts heat faster. This increase in temperature will have two effects: di- 
rect (hotter objects are brighter) and indirect (higher temperatures produce 
higher complex permittivities.) Therefore, for the multi-layer models, we 
calculate the surface temperature using a thermodynamic model based on 
certain assumed prevailing weather conditions. Since we assume the ice is 
in thermal equilibrium — valid if the ice is not too thick and the weather is 
changing relativly slowly — the temperature profile is linear. 
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Table 3: Thermodynamic model parameters for the two weather scenarios 





Scenario 


Parameter 


Fall freeze-up Winter cold snap 


Wind speed [m/s] 
Air temperature [K] 
Relative humidity 
Cloud-cover 
Insolation [W/m 2 ] 


2. 10. 
270. 260. 
0.5 0.1 
0.5 0.1 
50. 0. 



The following equation relates the ice surface temperature to net heat 



flux: 



hQ* =k(T w - T s 



0) 



where T w is the water temperature which is assumed to be constant at 
freezing (approximately -1.9° C at a water salinity of 35 psu), T s is surface 
temperature and k is the thermal conductivity of the ice. The net heat flux 
comprises the following components, with functional dependencies supplied: 



Q* = QeHT s )} + Q H (T S ) + QswiTf) + Qlw 



(10) 



The terms on the RHS are, from left to right: latent heat, sensible heat, 
longwave and shortwave; e(T) is the saturation vapour pressure. The first 
two terms are approximated with simple parameterisations while the long- 
wave flux is based on the Stefan-Boltzmann law. The shortwave flux is 
calculated primarily from geometric considerations based on the position of 
the Earth relative to the Sun. The following inputs are required for the 
model: surface- wind speed, humidity, air temperatu re and density (o r pres- 
sure), cloud cover and date and time or in solation. ( Cox and Weeksl . Il988l : 
Drucker et all . 120031 : lYu and Lindsavl . 120031 ) Inputs can be supplied to give a 
picture of the general weather conditions. For instance, fall freeze-up might 
be characterized by relatively mild temperatures, low winds, high humid- 
ity, high cloud cover and moderate insolation. By contrast, a winter cold 
snap would be characterized by low temperature, high winds, low humidity, 
clear conditions and little to no insolation. See Table [3l Equations ([9]) and 
(1101) a re solv e d wit h a numerical root-finding algorithm, specifically bisection 



Press et al.l . ll992h . 



The model can also be used as a crude ice growth model. The rate of 
ice growth is: 

where L is the latent heat of fusion for water. Empirical equations for 
determining the initial brine entrapment in sea ice have been derived by 
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Figure 8: Error residuals for the Vant effective permittivity models plotted 
against wavelength. Fitted curves are shown. 



Cox and Weeks! (|l988h and lNakawo and Sinhal (|l98lh and take the form 



S = So f{g) 



(12) 



where So is the salinity of the parent water and / is an empirical function 
of ice growth rate. 

As a qualitative validation of the salinity-thickness relations described in 
12.21 this model is quite effective. As seen in Figure [6j the model considerably 
overestimates the actual salinities, especially as the ice get thicker. This is 
because it does not take into account brine drainage and expulsion. 



2.4 Error 



Both the salinity-thi ckness relat i on fro m I Cox and Weeks! ()1974l ) and the 
mixture models from Vant et all ( 1978 ) are supplied with statistical error 
bounds (root-mean-square errors (RMSE) or residuals). We use these to 
generate error tolerances for the final 2Vthickness curves. Residuals from 
the Vant models were found to have a similar close correlation with wave- 
length as the relationship between the two effective permittivity models: see 
Figure This relationship was use to extrapol ate the residuals to higher 



frequencies. A Monte Carlo or "bootstrapping" (jPress et all 1 19921 ) method 
is used. Intitial salinities are perturbed with Gaussian deviates with zero 
mean and standard deviation equal to the residuals. So are the generated 
complex permittivities. The number of trials was chosen based on conver- 
gence of error estimates. 
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(a) Brightness temperature vs. ice thickness 
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b) Polarization difference vs. ice thickness 
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Figure 9: Modelled brightness temperature as a function of ice thickness 
using a single ice layer and constant temperature of 265K. Dotted lines 
show scattering component. 
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Figure 10: Comparis o n of modelled effective permittivities using 
Hoekstra and Cappillinol ( 197ll ) and the semi-empirical model described 



in the text. The Hoekstra and Cappillino mixture model was used by 
Naoki et al. (|2008h to model emissivity as a function of ice thickness. 
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(a) Brightness temperature vs. ice thickness 
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b) Polarization difference vs. ice thickness 




Figure 11: Modelled brightness temperature as a function of ice thickness 
using multiple (21) ice layers. Salinity profile has an 'S'-shape while the 
temperature profiles is determined with a thermodynamic model. Both are 
described in the text. Dotted lines show scattering component. 
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3 Results 



Here the final T^-thickness relations are presented. The error bounds are 
rigorous in the sense that those quantities with known error have been fully 
accounted for. However, there are many factors not accounted for in the 
model such as the possible effects of ice ridging and anisotropy of the brine 
inclusions. Factors without known error, such as the correlation length of 
the brine inclusions, have not been accounted for. 

The dependence on ice thickness of vertically polarized brightness tem- 
perature for the single (ice) layer model is shown in Figure along with 
error bars. Ice temperature was assumed to be 365 K. What these results 
show is a modest and fairly rapid increase in emissivity with thickness as 
the ice becomes progressively more opaque, thus less of the water shows 
through. Once the signal becomes "saturated" with only the ice emissions, 
the curve remains relatively flat, with an emissivity of close to one. The 
distance to saturat i on is controlled by the frequecy: from Equation (8) in 



Mills and Hevgsterl (|201ll ). the attenuation coefficient depends directly on 



the frequency. For all but the bottom two or three frequencies, the water 
layer hardly shows through even at the thinnest ice thickneses resulting in 
an almost constant emissivity. 

These results contrast with the those of Naoki et al. ( 20081 ) in which both 
modelled and measured ice emissivity was shown to have a strong and slow 
increase with thickness. The reason for the difference lies in the modelled 
effective permittivies, which are shown in Figure [TUl Naoki et al. use a mix- 
ture model that generates much larger values at higher salinities. It is this 
large real component that produces low emissities at thinner ice thicknesses. 
The changes in real permittivity will also produce higher polarisation differ- 
ences at th inner ice thickness and this pattern can be seen weakl y in Figure 
M(b). fSee lHevester et al.1 (120091 ) and lMills and Hevgsterl (|201lh for an ex- 
planation.) It is this decrea se in polariza t ion d ifference that is frequently 



used to detect thin ice, as in iMartin et all (|2005h 



The 2Vthickness curves for the more sophisticated, multi-layer model are 
shown in Figure [TT1 21 layers were used in the model. Note that this model 
shows a stronger increase/decrease in brightness-temperature/polarization 
difference with ice thickness, especially at higher frequencies. This is due 
mainly to the scattering component, which is shown by the dotted lines. 
Error bars for these curves are narrower because the random variations in 
effective permittivity in each layer will tend to cancel. Note that errors 
are higher for the polarization differences (in both models) because they 
were calculated from the root-sum-square of the errors in the horizontal and 
vertical. 
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4 Summary and conclusions 



Radiative transfer-based emissivity models were run for all frequencies of the 
Advanced Microwave Scanning Radiometer on EOS (AMSR-E) instrument 
(6.925, 10.65, 18.7, 23.8, 36.5 and 89.0 GHz) as well as 1.4 GHz in order to 
determine the relationship of emissivity to sea ice thickness. Because of ice 
growth processes, ice s alinity is normal l y inv ersely related to ice thickness. 
The satistical model of Cox and Weeks ( 1974! ) was used to relate bulk salin- 
ity to ice thickness. The emissivity model assumes plane parallel geometry; 
two versions were employed: one with only one ice layer (uniform proper- 
ties throughout) and one with multiple layers. For the single layer model, 
a constant temperature of 265 K was used. For the multi-l ayer model, th e 
salinity profile was modelled with the 'S'-shaped profile from Eicken ( 19921 ). 
scaled to match the bulk salinity. For the temperature profile, a thermody- 
namic model was employed which assumed thermal equilibrium and weather 
conditions equivalent to fall freeze- up. 

Complex effective permittivitie s to feed to t he RT model were calculated 
from salinities and temperatures. Vant et al. ( 19781 ) derived measurement- 
based, statistical models that linearly relate effective permittivity to brine 
volume. These models were only generated for frequencies between 0.1 and 
4 GHz. In order to extend t he Vant models to high e r freq uencies, theoretical 
mixture models taken from Sihvola and au Kong ( 19881 ) and based on the 
low frequency limit were statistically adjusted to the Vant models. 

Error tolerances were calculated using a Monte Carlo method based on 
the residuals supplied with both the Cox and Weeks salinity-thickness rela- 
tion and with the Vant effective permittivity models. 

When scattering is neglected, the brightness temperature is found to vary 
only slightly with ice thickness for both the single- and multi-layer models, 
except at the lower frequecies, where brightness temperature increases with 
ice thickness. The latter effect is caused by the greater translucency of the 
ice at lower frequencies: the water below shows through. The reason little 
change is observed at higher frequencies is because the effective permittivity 
models show only a we ak increase i n real permittivity with brine volume. 
Contrast this result to Naoki et al. ( 20081 ). Note that both the real and 



imaginary permittivities are almost linearly related to salinity. 

Adding scattering produces much more variation with ice thickness of 
both brightness temperature and polarization difference, particularly in the 
multi-layer model. Below about 30 cm, the brightness temperature is seen 
to increase sharply with ice thickness, while above that it decreases gently, 
especially at the higher frequencies. This is in accordance with measure- 
ments where it is observed that thin ice has both a lower brightness tem- 
perature and higher polarization difference, while older, thicker ice tends to 
be radiometrically cooler. Meanwhile, the size of the scatterers ensure that 
scattering is greater at higher frequencies. It also makes sense that most of 
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this difference is due to scattering: thinner ice is more saline, hence it will 
have both more and larger brine pockets. It also tends to be more granular, 
although this is not well accounted for in the models. For the multi-layer 
model, scattering in the horizontal polarization is larger than that in the 
vertical, a possibly questionable result. 

One puzzling feature of the muliti-layer results is the low polarization 
difference for very thin ice (below 20 cm), especially at higher frequencies. 
This is due to the scattering component and does not accord well with 
measurements. To correct this, a better understanding of the relationship 
between scattering and ice physical properties is necessary. 

Error bounds were larger for the single-layer model and reached as high 
as 50 K. For the multi-layer model, bounds high for thin ice thicknesses 
(again, as high as 50 K) but tended to be lower. Most were below 10 K, 
particularly for lower frequencies and higher ice thicknesses. 

Better ice emissivity models require better estimates of ice effective per- 
mittivity, hence the design of a new, combined model. The correlation in 
wavelength of the Vant empirical models with the Sihvola and Kong the- 
oretical models (based on the low- frequency limit) is likely significant, but 
a fuller understanding would require a deeper investment into the theory. 
A similar statement is true for the scattering components of the model in 
which the correlation length was related to brine volume using a simple, ad 
hoc assumption. 

Understanding ice growth processes would be another area for future 
study since they are closely connected with ice emissivity. A simple ice 
growth model based on the thermodynamic model for ice temperature was 
used to qualitatively confirm the relationship of bulk salinity to thickness. 
Ice cores from the Weddell Sea were similarly modelled however the results 
are not shown here. 
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